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Abstract 

A mean-field formulation is used to investigate the bifurcation diagram for 
globally coupled tent maps by means of an analytical approach. It is shown 
that the period doubling sequence of the single site map induces a continuous 
family of periodic states in the coupled system. This type of collective motion 
breaks the ergodicity of the coupled map lattice. The stability analysis suggests 
that these states are stable for weak coupling strength but opens the possibility 
for more complicated types of motion in the regime of moderate coupling. 
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1 Introduction 

Despite the deep insight which has been gained into the dynamics of low dimensional 
chaotic systems and the sophisticated methods that have been developed for its 
analysis (cf.|l|, there is a tremendous lack in the understanding of the influence of 
irregular motion on high dimensional dynamics. The failure is twofold. On the one 
hand it is not quite clear what kind of quantities are suitable for the characterization 
of chaotic motion in extended systems. On the other hand there is no common sense 
what kind of model systems capture the relevant features of chaotic motion in high 
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dimensional phase spaces. Concerning the second problem systems of coupled maps 
have been proposed as suitable models |§ because of the simplicity of numerical 
simulations and of the fact that maps have proven as reasonable models for the 
investigation of low dimensional chaos. Although there does not exists a satisfactory 
derivation of coupled map lattices form physically sound equations of motion a rough 
estimate suggests Q that they capture important features of realistic systems[]. In 
addition to a pure numerical simulation, analytical and even rigorous approaches 
have been applied successfully to simple coupled map lattices. Especially it has been 
shown that for rapidly decreasing interaction between the different lattices sites the 
hyperbolic property of the local dynamics is inherited to the coupled system, so that 
the space-time correlations decrease exponentially || |J. This mixing property is 
sometimes used as an operating definition for space-time chaos although it does not 
seem to be fully satisfactory. It has been suggested that a breakdown of this regime 
can be understood in terms of phase transitions in the two dimensional spin lattice 
which serves as a symbolic dynamics 0. 

In order to contribute by means of an analytical approach to the problem how 
chaos influences the motion in high dimensional phase spaces I will consider in this 
article a rather simple but nontrival coupled map lattice. For that purpose tent maps 
are chosen, as they have nontrivial features concerning their bifurcation scenario, 
their symbolic dynamics and their spectral properties but allow for an analytical 
treatment (cf. ||). For the coupling mechanism an all to all interaction is the simplest 
choice so that one ends up with 

(i-£)/(4"») + £ A„=:r„(4")) (i) 

/j=0 

1 — a\x\ . (3) 

Here N denotes the system size and we will dwell on the limit of large system size 
in the sequel. A superficial inspection of the central limit theorem together with the 
chaotic properties of the local dynamics would suggest that global quantities like the 
coupling field h n possess fluctuations which decrease as iV -1 / 2 with the system size. 
But such a property requires the factorization of the spatial correlation function at 
different sites (cf. M) and is in this sense equivalent to the mixing property referred 
above. It may be violated due to the infinite range of coupling^, so that global 
quantities obey a nontrivial dynamics even in the limit of infinite system size ||10|| . It 
has been suggested that the occurence of such a kind of motion, which corresponds 
to a partial synchronization in the extended system, is related to the nonhyperbolic 

1 Contrary to some suggestions which can be found in the literature the coupling does not derive 
from a discrete Laplacian but from the propagator of the extended system. 

2 Sometimes the slightly misleading expression "Violation of the law of large numbers" is used 
for this phenomenon in the literature. 
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properties of the local dynamics [JTTJ] . For the simple model system ([T])-([j|) I want to 
investigate how such a kind of motion occurs. 

To keep the presentation 



For that purpose a mean-field like approach is used [12 



self-contained and to set some notations this formulation will be recapitulated in 
section 2. With its help the bifurcation behaviour of periodic densities is analysed in 
section 3. The highly nontrivial aspect of the stability of these solutions is addressed 
in section 4. 



2 Mean— field formulation 

The value of the coupling field (0) is solely determined by the "one-site" density 
by dint of the relation 



N 



f(x)p n (x)dx . (5) 

If one recalls that the dynamics is governed by the effective map T n then the time 
evolution of the density (f|) is given by 

Pn+\{x) = J 5(x- T n (y))p n (y) dy =: C n p n {x) (6) 

where the map T n depends itself on the density (cf. eqs.(|J), @). In this respect 
eq.(HP constitutes a nonlinear evolution equation for the "one-site" density It 
should be stressed that this expression is an exact consequence of the full evolution 
equation ([]]) as long as the density takes the form (|]). In additon the system size 
does not enter the formulation in an explicit way. 

It is supposed that the limit of infinite system size can be performed in such a way 
that instead of the special form (|j) every integrable density is accessible during the 
evolution. Although this limit, understood in the weak sense, seems to be difficult 
to justify on a rigorous basis this assumption is highly plausible on physical grounds 
and can be supported by numerical simulations. By this assumption those transients 
which increase with the system size are turned into stable solutions of the mean- 
field description^. Hence the stationary solutions of eq.(^|) capture also the transient 
behaviour which becomes quasi-stationary in the limit of large system size and which 
governs the dynamics of large systems on accessible time scales. 

The density (§) allows for the computation of any globally averaged quantity. In 



addition it has been demonstrated ||T3| that in cases where the the mean-field dy- 
namics (§) converges to a stable solution, the correlation properties between different 
sites can be obtained from the "one-site" density also. 

3 It is a simple task to demonstrate this property on the example of globally coupled shift maps. 
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3 Periodic Densities 



We are going to construct periodic solutions po, . . . , p p -i, p p = Po of the evolution 
equation (^) by generalizing some ideas of jnj . For reasons that will become obvious 
in the sequel we consider periods p = 2 M . Eq.([|) tells us that p is determined by 
an invariant density of the map 

(T p _i o . . . o T ) (x) = r? p _i - a red \r] p _2 - ... - a red \r] - a red \x\ | . . . | . (7) 

Here the abbreviations 

a red := (1 - e)a, r) n := 1 - ea J \x\p n (x) dx (8) 

have been introduced. rj n is positive as long as the coupling strength satisfies e < 1. 
We further note the trivial fact that every mean-field map Ti(x) = T]i — a red \x\ is 
topological equivalent to a single tent map with a reduced parameter 

F{z) = 1 -a red \z\ (9) 

by a linear coordinate transformation x = T]iZ. 

Let us first consider the simple case that all parameters (H) are equal i] = . . . = 
Vp = V*- The p iterate of the map (|9|) admits of p different ergodic components 
for V2 < a red < v2. They correspond to the p-band chaos in the tent map. Let 
Vi(x),i = 0, . . . ,p — 1 denote the different normalized ergodic absolutely continuous 
distributions and let uq denote that ergodic component which contains the critical 
point z — 0. They obey the Frobenius-Perron equation which reads Vi + i(x/r]*) = 
£*^j(x/?7*). For l^z^p— lit just leads to an affine transformation of the 
distributions v^. The density 

Mx) ' (10) 

constitutes by construction a p periodic orbit of eq.@ which in fact is a fixed point. 
Finally the free parameter 77* is determined by the relation (R) 



V* = 1 ~ 7^1 / dz (11) 

where it is worth to mention again that the distributions Ui depend only on a red by 
definition. The meaning of this fixed point solution is quite simple. Every periodic 
chaotic band contributes to the density fllOD in such a way that the time dependencies 
cancel each other exactly. Finally let us denote for later reference by <jj the symbol 
sequence of the critical point that means <7j = sgn(T^ +1 (0)). 

We continue our construction of a continuous family of periodic solutions by 
giving the different ergodic components different weight factors. For that purpose 
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Figure 1: Diagrammatic view of the second iterate T\ o Tq for \[2 < a re d < \/2 and two 
different choices of parameters rji. 



one has to show that the components are not destroyed if the parameters r] n are not 
equal. In fact as sketched in Fig. [l| the change of these quantities will lead only to 
an affine transformation of the different components but not to topological changes. 

To make this statement explicit consider a p periodic sequence of parameters rji 
close to the value r/* mentioned above. Fix a neighbourhood of the critical point 
x = that contains the support of the former ergodic component. As the finite 
itinerary of the critical point a , . . . , a p is unaltered for rji sufficiently close to 77* the 
iterated map (0) on this neighbourhood reads 

p-1 p-2 p-2 

(T p _ 1 o...oT )(x)=^(-a re ,) p -'- 1 -]^a fc -r /z + (-a red ) p - 1 -^[a fc -|x| . (12) 

1=0 k=l k=0 

Since this expression is in a neighbourhood of the critical point topological equivalent 
to the p th iterate of the single tent map (|D, eq.(|l2|) admits of an ergodic component 
^o(^)- It is given in terms of the corresponding ergodic component vq of the single 
map by 

z> (x) = -v (pj (13) 
where the scaling factor T follows from the linear conjugacy to beQ 

p— 1 p— 2 p— 1 p—2 

r^-twr'- 1 ■x\vk = Y;(- a redT 1 - 1 -H^k-vi ■ (14) 

1=0 k=l 1=0 k=l 

4 As long as the itinerant is not periodic the sum on the left hand side does not vanish. 
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By construction the sequence of densities := LiVi is periodic with period p. In 
view of the fact that the mean-field map Tj is conjugate to the single tent map @, 
and that the former acts as an affine transformation for rji sufficiently close to 77*, it 
is a simple task to express the densites z>j in terms of the ergodic components z/j of 
the single map 

= ^i(^) (15) 



r V r 

i-1 i-2 

5 t := ^(-a^-^-JJ^'^-r) • (16) 

1=0 k=l 

Eq.(|15D just states that the densities follow by an affine transformation from the 
ergodic components z/j. In order that this representation is valid it is necessary, that 
the support of these densities is not shifted across the critical point, that means <5j has 
to be sufficiently small. But from eqs.flHD and ( ]T5| ) it is obvious that this property 
holds for r/i sufficiently close to 77*. Now by construction the p periodic convex sum 

p— 1 p— 1 

Pn{x) = 22a n +ih(x), a n+p = a n e [0, 1], = 1 (17) 

i=0 k=0 

satisfies eq.@. It remains to show that the weights ctj can be chosen in such a way 
that the self-consistency condition (|8|) is valid. For that purpose one inserts the 
density (|17D into eq.(|[) and obtains 

P-i 

7] n = 1 - ea a n+i (TQ + (7i-iSi) . (18) 

i=0 

Here we have expressed the integrals in terms of the densities z/j of the single map 
and have used the abbreviation 



Ct := J \z\ui(z)dz . (19) 

But now eqs.flHf) and the definitions (0) and ( |IB| ) yield a system of p linear equations 
for the p coefficients 7] n if one considers the weights as given. The solution of this 
system completes the construction of the periodic density^. It is completely given in 
terms of the densities of the single map (|9|). 

Hence one has found a continuous family of p = 2 M periodic solutions for y2 < 
a re d < ^2 which emerges by attributing slightly different weights to the different 
ergodic components of the iterated mean-field map (cf. eq.([l7])). Fig. |2| summarizes 
this result in a partial bifurcation diagram. The period doubling of chaotic attractors 

5 One might argue, that the system is singular for a large set of parameter values (a, e). But as 
for fixed value of a re d the coefficients depend on the parameter ea in a linear way the determinant 
vanishes at most at p isolated points. 
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Figure 2: Bifurcation diagram for globally coupled tent maps. Pn denotes the region 
where a continuous family of n periodic densities exists. The shaded area indicates the 
region where asymptotic periodicity gets lost. 



of the single map translates into the bifurcation of periodic densities of the coupled 
system. At each bifurcation line the period p family of densities splits into a period 
2p family of densities. As was already mentioned in [|ll| at a re d = 1 asymptotic 
periodicty gets lost. In addition the bifurcation at a re d = \/2 has been observed in 
numerical simulations some time ago |L2] . 

I have yet not claimed that the periodic densities are stable. One might argue 
that on physical grouds these solutions are stable, as the ergodic components are 
preserved in the coupled map lattice. But as was already mentioned in |nj the 
problem of stability turns out to be highly nontrivial. Furthermore in numerical 
simulations of finite lattices one observes additonal oscillations ITBi which cannot be 



attributed to the period doubling scenario. Before I turn to the discussion of the 
stability let me demonstrate by means of a numerical approach that the densities 
might be stable at least for small coupling. For that purpose simulations of eqs. ([[])- 
(§) have been performed for parameter values within the period 2 and 4 window. 
For different lattice sizes one finds a corresponding periodic motion in the coupling 
field on which a "noisy" component is superimposed. In order to check whether this 
component is a finite size effect the mean square deviation, that means the second 
cumulant of the p th iterate of h n , has been analysed in dependence on the system 
size (cf. Fig.^|). The behaviour according to the central limit theorem suggests that 
the nonperiodic component is a finite size effect and that these periodic states are 
dynamically stable. But the numerical values of the periodic density depend on the 
initial condition. Hence the latter determines which member of the continous family 
of periodic densities is selected as the stationary state. 
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Figure 3: Mean square deviation A of the second respectively fourth iterate of the coupling 
field h n for e = 0.1 and two values of a re( i in dependence on the system size. 



4 Stability Analysis 

The linear stability analysis of the equation of motion (|6]) faces the problem that 
the densities are in general not smooth functions. Following the ideas of JDJ this 
equation takes a somewhat simpler form if one restricts the analysis to the case of 
densities which are made of a countable infinite number of step functions 



3=0 



(20) 



Due to the piecewise linear shape of the tent map this form is preserved during the 
evolution. Then eqs.@, (||) read (cf. [ fPT| ) 



b° 

°n+l 



°n+l 



c n+l 

T n {x) 



3=0 



1 



-^sgn(4) 

^red 

T n (0) 

1 - e + eh n - a red \x\ 

3=0 



(21) 

(22) 

(23) 
(24) 
(25) 

(26) 



We will consider only the stability properties of the fixed point solution fllCf ) 
in the sequel. In terms of the representation (pH) this solution is obtained from 
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eqs.(0)-(]2|) 



ared/ k=0 

ci = Tt\0) (28) 

jr=0 V rea/ fc=0 

where again crj := sgn(e£) denotes the symbol sequence of the critical point. To 
investigate the stability properties of this solution let us consider the eigenvalue 
problem that emerges from the formal linearization of the evolution equations (|2"T|)- 



A/5° = — + (30) 

3=0 

= —^(Tj (31) 

A 7 ° = e5h (32) 
A7 J+1 = e5h — a r& pf 3 Oj (33) 

5/1 := "fEWl^ j + 2 VM\l j ) ■ (34) 

3=0 

Since these equations decouple the eigenvalue problem can be analysed on both 
invariant subspaces separately 
With the solution of eq.([3l|) 

^-^)>1>* < 35 > 

\ rea / k=Q 

eq.flSTf) yields the characteristic equation 

1 

= (1-z) J] z:=— . (36) 



j=0 fe=0 



The sum on the right hand side is just the A-expansion and yields the spectrum of 
the Frobenius-Perron operator of the fixed point map T* |§. Hence the spectrum is 
contained within the unit circle and the eigendirections belong to the center-stable 
manifold of the fixed point density Obviously there occurs a doubly degenerated 
eigenvalue A = 1. It is caused by two constants of motion . W n = 0, J2j Hi c n = 1 
of eqs. (]26|) which reflect the normalization of the density and the boundary 
condition p n (— 1) = 0. 
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Using the solution of eqs.(p2|), (|33| 



^ = Eri(-T*"^)-x (37) 

1=0 k=l 

eq.(Q) yields for ft = the characteristic equation 

i = -y£(-^)'£ft:k:ilp- 1 . m 

1=0 j=l k=l 

If one takes the representation ([27]) of the fixed point into account one ends up with 

z=o 



where the coefficients are given by 



* :=&££(-— Vl^lff** • (40) 

J=0 x 7 fc=0 



Although one has obtained a quite simple expression which can be evaluated by 
means of a numerical approach^ it seems to be difficult to make general statements 
about the solutions of eq.([40"D respectively its analytical continuation. But if one 
restricts to the simple case that the itinerary terminates in the fixed point of the 
map T*, that means \c[\ = Cp, I ^ l , then the coefficients obey 

ai = b°c F J2(-—) J li° k = ^ l>h ■ (41) 

3=0 \ ared ' k=0 

To establish this identity one has used the fact that the sum represents the A- 
expansion of the fixed point map T* and that its Frobenius-Perron operator admits 
of the eigenvalue A = 1 |§. Hence the characteristic equation fl39|) reduces to a 
polynomial which yields for sufficiently small coupling only eigenvalues within the 
unit circle. 



5 Discussion 

The bifurcation diagram depicted in Fig. |2] summarizes main results of this article. 
It provides a foliation of the region of asymptotic periodicity a re d > 1. In each 
region a continuous family of periodic states, which is generated by the different 



3 The crucial step in this approach is the desired truncation of the potentially infinite series. 
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ergodic components of the iterated tent map, accounts for the time dependency of 
the stationary solution. 

The stability analysis suggests that these states are stable at least for small cou- 
pling strength e. But there are several limitations in this analysis. For the stability 
of the fixed point density the estimates presented in section 4 are not uniformly valid 
in the parameter space and the problem whether stability is attained on a parameter 
set of sufficiently large measure is not solved. However numerical simulations and 
a numerical analysis of the eigenvalue equations (|30|)- (|34T ) [ [Dfl suggest stabilty for 
small coupling. In addition one should recall that the discussion of stability pre- 
sented above was based on a formal linearization of the evolution equation. Such a 
procedure takes deviations of the parameters c J with equal weight into account. This 
approach seems to overestimate the effect of these deviations on the density fl20|) 
since the coefficients V decay exponentially. Hence one has to specify the space on 
which the eigenvalue problem has to be treated which may restrict the spectrum to 
a few relevant eigenvalues (cf. e.g. [0 for such a phenomenon in the context of the 
Frobenius-Perron operator). Despite these shortcommings the approach presented 
here might be a basis for further investigations. 

The periodic states investigated in this article explain in a clear fashion the time 
dependence of globally averaged quantities. Therefore one mechanism for time fluc- 
tuations in globally averaged quantities seems to be well understood. But numerical 
simulations of coupled tent maps show additional quasiperiodic solutions fIH[] . The 
explanation of this kind of collective motion which resembles in some respect a Hopf 
bifurcation seems to be one challenge for the further research in this field. In addi- 
tion it might be this kind of motion which appears in the frequently analysed case 
of coupled logistic maps. 

The feasibility of analytical computations on piecewise linear map lattices make 
these kind of models suitable for studying the influence of chaotic motion in high 
dimensional systems. Even the apparently simple model treated in this article shows 
several unexpected features and it is far from being solved completely. Other prob- 
lems may be attacked by this kind of model also, e.g. the effect of a long but not 
infinte range of coupling in connection with the validity of a mean-field treatment, 
or the construction of the symbolic dynamics and the effect of pruning of symbol 
sequences. But these questions are left for future work. 
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